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Abstract 

In this short paper, the validity of Koopmans' theorem in the 
Hartree-Fock theory at non-zero temperature (Hartree-Fock statistical 
theory) is investigated. It is shown that Koopmans' theorem does not 
apply in the grand-canonical ensemble, due to a missing contribution 
to the energy proportional to the interaction between two electrons 
belonging to the same orbital. Hartree-Fock statistical theory has also 
been applied in the canonical ensemble [Blenski et al., Phys. Rev. E 
55, R4889 (1997)] for the purpose of photo-absorption calculations. In 
that case, the Hartree-Fock self-consistcnt-field equations are derived 
in the super-configuration approximation. It is shown that Koopmans' 
theorem does not hold in the canonical ensemble, but that a restricted 
version of the theorem can be obtained, by assuming that a particular 
quantity multiplying the interaction matrix element in the expression 
of the energy does not change during the removal of an electron. 



1 Introduction 

In atomic physics of hot dense plasmas, the Hartree-Fock (HF) theory is 
applied to model self-consistent potentials of atoms and ions. Another in- 
terest of HF theory is to use HF one-electron states in order to calculate the 
average configuration energies. As shown by Mermin [I], the HF theory at 
non-zero temperature is equivalent to minimizing the grand thermodynamic 
potential in a restricted class of statistical operators. In the following, such 
a theory will be refered to as "grand-canonical statistical HF theory" . 

A few years ago, Blenski et al. [21 [3] generalized the statistical Hartree- 
Fock approach to atoms in plasmas at finite temperature in the framework 
of the super-configuration approximation successfully applied in the STA 
(Super Transition Array) method [4]. In the following, such a theory will be 
refered to as "canonical statistical HF theory". 
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Koopmans' theorem [21E1I7] states that, when the number of electrons qi 
in a given orbital i is greater than 1, the interaction within the group is gen- 
erally much smaller than the energy difference resulting from the removal of 
an electron. Therefore, e» represents roughly the average theoretical ioniza- 
tion potential for a single electron during the removal of the group. Koop- 
mans [5] identified the physically meaningful solutions as those for which 
was still an approximation to the ionization potential when = 1, and an 
approximation to the average ionization potential otherwise. An important 
point is that Koopmans' theorem is valid only if the one-electron wavefunc- 
tions in the iV-electron and the (iV±l)-electron Slater determinants are the 
same. 

The purpose of this short paper is to check whether Koopmans' theorem 
is valid in the framework of statistical HF theory in the grand-canonical 
ensemble [JJ and in the canonical ensemble [2j [3]. For that purpose, the 
energy of the system with qi electrons in orbital i is compared to the one 
with qi — 1 electrons in orbital i. The difference is related to the value of 
the energy of orbital i. 

2 Grand-canonical ensemble 

2.1 Average energies 

The total energy of a configuration with N orbitals (raiZi), (712/2), ■ • • , (i^nIn) 
is i 

N N 
E = J2qi + qi {qj - S itj ) V tj (1) 

i=l i,j=l 

where qi is the population (number of electrons) of orbital i and, in the 
non-relativistic case (Schrodinger equation), the one-electron energy Ii reads 
(atomic units are used throughout the article): 

where Pi(r) is the radial part of the wavefunction multiplied by r. The 
interaction matrix elements are 

and for % ^ j, 
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where and are the direct and exchange Slater integrals respec- 
tively [Qj. The quantity ^ m re P resen ^ s a 3j-symbol [10J. The 
average energy of the configuration is 



-. 91 92 93 9N 



N , AT 



2 

i=l *J=1 



91=0 g2=0 g3=0 9jv=0 
AT / , 

x II ( ^ ) e-^-ri Qk (5) 
k=i \ qk ' 

where \^\ = a '/(^K a ~~ &)0 i s the binomial coefficient, /i the chemical 
potential and Z the partition function of the system: 

91 92 93 9iV JV / \ 

gi=0g 2 =0 93 =0 gjv=0fc=l V Hk J 

One finds that the average energy flS} can be written 



iV N 



i=i j,j=i ^ ^ 

AT AT 

= ^2(qi) J i + o («) (to) ~ Pi Vi i> ( 7 ) 



2 

i=i ij=i 



with 



(gi)=giPi= 1 + e _ g ; (ei _ M) - (8) 

The energy required in order to remove one electron from orbital i is 
defined by 

A£« = E ({q k - 5 ik , k = 1, AT}) - fc = 1, N}) , (9) 

where <5j& denotes Kronecker's symbol. The notation {q k — 5ik,k = 1, iV} 
means that the population of orbital i is reduced by 1, the other populations 
remaining unchanged. Therefore, one obtains 

A£ (i) = -h - (<«> - ^ (l - ±) Vfc - (10) 
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2.2 Expression of orbital energies 

The Hartree-Fock equations read [6] 





" d? i(i + iy 


H 


dr 2 r 2 



^ + V{r) + 5Vi{r)-ei J Pad +„,{,■) a (1.1.) 



with 

N 



J ' poo i 

V(0=E<W ^(rO-PiCr') dr'. (12) 



The correction <5V^ is given by 



SV i (r) = -PiV^(i,i)-^^J2( l o o S)V(*)(i<) ((*>-«), (13) 



where 



and 



= [°° Pj(r')^Pi(r') dr' (14) 
Jo /> 



^) = -E¥E(q S ^^^OPiW. (15) 
Multiplying equation (jlip by Pj (r) and integrating from to oo leads to 

e l = h + («) f 1 - -) V u + E(*)^i- (16) 
Therefore, one has 

A£« = - €f + \ (l - -) Vu ± -e h (17) 
2 V 9iJ 

which means that Koopmans' theorem is not verified in that case, be- 
cause of the remaining term (1 — l/gi)Vu/2. Each orbital is shifted from 
that quantity, which is the signature of the competition between the di- 
agonal term of the interaction matrix and the degeneracy. We can see in 
table [H in the case of a carbon plasma at T=30 eV and p=0.01 g/cm 3 , that 
the difference between the energy difference AE^ and e« can reach 25 %. 
The shift is more pronounced for the higher-energy orbitals, which can be 
explained by the fact that the electrons in such orbitals are very sensitive 
to electron-electron interactions, unlike electrons in the lower orbitals which 
are more subject to the attraction of the nucleus^. 



2 However, it is worth mentioning that Janak's theorem [111 118] holds in that case, i.e. 
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Orbital 


9i 


Energy t\ (eV) 


V u 


(1 - l/g l )V ii /2 


AEW (eV) 


Shift (%) 


Is 


2 


-383.376 


44.892 


11.223 


394.599 


2.93 


2s 


2 


-57.101 


9.212 


2.303 


59.404 


4.03 


2p 


6 


-48.905 


15.994 


6.664 


55.569 


13.63 


3s 


2 


-17.519 


3.791 


0.948 


18.467 


5.41 


3p 


6 


-16.996 


6.312 


2.630 


19.626 


15.47 


3d 


10 


-16.188 


8.088 


3.640 


19.828 


22.48 


4s 


2 


-6.934 


2.062 


0.515 


7.449 


7.43 


4p 


6 


-6.154 


3.749 


1.562 


7.716 


25.38 



Table 1: Energies (in eV) of orbitals Is to 4p for a carbon plasma at T=30 
eV and p=0.01 g/cm 3 . The chemical potential is ^=-186.702 eV. 



3 Canonical ensemble 

The formalism presented in the preceeding section allows non-integer pop- 
ulations for the different orbitals. Therefore, it provides the average (in the 
sense of the most probable) configuration of the plasma. From this average 
configuration, the real configurations can be built, by rounding the popu- 
lation values to the closest integer. The problem is that, in hot plasmas, 
the number of configurations can be really huge. The super-configuration 
method [I] has been invented in order to remedy this problem. A con- 
figuration is made of orbitals with integer populations; in the same way, 
a super-configuration is made of super-orbitals with integer populations, 
a super-orbital being a group of orbitals which energies are close to each 
other. For instance, (Is2s2p) 3 (3s3p) 4 (3ci) 7 is a super-configuration made of 
3 super-orbitals populated respectively with 3, 4 and 7 electrons. 

The super-configuration approximation enables one to calculate the equa- 
tion of state beyond the average-atom model [121 E2 El- The idea is to 
study the influence of the population fluctuations on the thermodynamic 
quantities. For instance, the pressure of the plasma is given by 

P = ]TW H P H , (19) 

where W= and P= represent respectively the probability and the pressure 
of super-configuration 3. In the particular case where a super-configuration 
is an ordinary configuration, the equations reduce to the standard HF equa- 
tions. The finite-temperature HF method was derived in the framework of 



d(cn) 
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the super-configuration approximation by Blenski et al. [2J [3] . The aver- 
age shell populations and interaction matrices, which are averages of the 
corresponding quantities for configurations, are given in terms of partition 
functions [15} \W[ I17j . This allows one to avoid problems stemming from 
non-integer population numbers in other thermal HF theories |18j . 

3.1 Expression of the average energy of a Q-electron config- 
uration 

The Hartree-Fock equations read [3] : 
d 2 1(1 + 1) 



dr 2 r 



2 



^ + V(r) + 5V i (r)-e i j P,(r) +//,(/•) 



1 

N 

= ^ S k,h P ^)- ( 2 °) 

The different terms in the left-hand-side of Eq. pip are 



JV 

5Vi(r) = ^ {5 itS (g { - 1) (Si - Pi) 5 kfi - Si] 

s,k=l 

+ (1 - S itS ) (q s ) Q (g s - 1) H is 6k,o}Y$(r), (21) 

and 

N N 

hi(r) = - Yl (q s ) Q [H is + l]Ps(r)J2 Y sf(r), (22) 

s=l,s^i k=l 

with 



.(h), > _ 9i Ik k l s " 



Y " (r) = ^S~) { ) /- ^^W) *>. (23) 



o rv 



There is a small typographical error in Ref. [2j: a factor (g s — 1) is 
missing in Eq. (6c). However, in Eq. (7c) of Ref. [3], the expression is 
correct. Multiplying equation (|20p by Pj(r) and integrating over r from 
to infinity leads to the orbital energy 

€i = Ii + Si ( 9i - 1) Vu + ^ [1 + H^] Vis(q s ) Q . (24) 

Let us consider the following super-configuration with N orbitals and Q 
electrons: 
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(mh n 2 h ■ ■ ■ n N l N ) Q . 
Its average energy reads 

N j N 

2 

i=i 

where 



(?i>Q = 9iPi = ~9i (~ Xi 



71=1 



UQ-n(g) 

U Q (g) 



is the average population of orbital i, = e ^ ^, and 

Si 92 gz 9n N 

k 



31=0 c?2=0 i}3=0 <?jv=0fc=l 



One has also 

W rs = 1 + (5 rjS 



EAT „ 
j=l*= 



>gr (gr ~ 1) 
(9r)Q - 1 



- 1 



where 



1 Q 

Si = -E( n - 1 )(- Jf < 

W n=l 



+ (1 - S r ,s) H rs , 



,n U Q-n(g) 



Uq{ 9 ) 



and 



X s — X r 



5s -^s <?r -X'r 



(9 S }Q (9r>Q 



- 1. 



(25) 



TV j iV 

(^)Q = J>i>Q J '+oE ^>Q (^)q - Pi 5m) Wij (26) 



(27) 



(28) 



(29) 



(30) 



(31) 



The energy required in order to remove one electron from orbital i is 



A(£«) Q = E ({(q k ) Q - 5 ik , k = 1, N}) - E ({(q k ) Q , k = 1, N}) + -e u 

(32) 

which means that Koopmans' theorem does not hold in the statistical 
Hartree-Fock theory for the canonical thermodynamic description of the 
system. If Wij did not change during the removal of an electron belonging 
to orbital i, one would have 

A£«- = I i -S t (g t - 1) V u -£[1 + H is\ Vis(q s )Q = -<*, (33) 

but this is not true in general. 
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4 Conclusion 



Koopmans' theorem does not apply in the statistical Hartree-Fock theory 
neither in the canonial ensemble, nor in the grand-canonical ensemble. In the 
grand-canonical ensemble, an additional term exists in the energy variation 
due to the removal of an electron from orbital i. It represents a shift of the 
orbital energy depending on its degeneracy and on the diagonal matrix ele- 
ment describing the interaction between two electrons in that orbital. In the 
canonical ensemble, one finds, in the framework of the super-configuration 
method, that Koopmans' theorem does not hold stricto sensu, unless a par- 
ticular quantity is ensured to be unchanged when the number of electrons 
is decreased by one. 
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